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INTRODUCTION 



HYMO (7) 3 is a problem-oriented computer 
language for modeling surface runoff and sedi- 
ment yield from watersheds. The language is 
called HYMO from the words "hydrologic 
model." HYMO was designed for planning flood 
prevention projects, forecasting floods, and re- 
search studies. It consists of a main program and 
16 subroutines written in FORTRAN IV, but it 
can be used by hydrologists with little knowl- 
edge of computer programing. The language 
provides 17 commands for the hydrologist to 
use in any sequence for application to any 
watershed. 

HYMO was designed to transform rainfall 
data into runoff hydrographs and to route these 
hydrographs through streams and valleys or 
reservoirs. It will also compute the amount of 
sediment produced by a storm at any point on a 
watershed. It will be useful to research hydrolo- 
gists in studying the effects of watershed and 
storm characteristics on the flood hydrograph. 
HYMO is also a good research tool for testing 



hydrologic procedures; for example, a new 
flood-routing method could be added to HYMO 
and tested easily, because the inflow hydro- 
graphs and the rating curves are available in a 
HYMO program. 

HYMO is flexible. Present hydrologic proce- 
dures can be modified or deleted, and other 
hydrologic procedures can be added by hydrolo- 
gists familiar with FORTRAN IV programing. 
Adding a new command simply requires the 
addition of a new subroutine. 

HYMO is efficient, practical, and generally 
applicable. HYMO programs can be written and 
the results interpreted by hydrologists who have 
no conventional computer programing experi- 
ence. The hydrologic procedures used in HYMO 
are practical required inputs are easy to obtain 
for most watersheds. 

HYMO was written for the IBM 360-65 
computer, but it could be run on an IBM 1130 
with little modification. The storage require- 
ment is about 73 K. 



OPERATION OF HYMO 



HYMO consists of a main program and 16 
subroutines. The HYMO card deck is set up in 
the following order: 



'Hydraulic engineer, College Station Area, Southern Region, 
Agricultural Research Service, U. S. Department of Agri- 
culture, Riesel, Tex. 

= Head, Environmental Engineering Division, Civil Engineer- 
ing, Texas A&M University, College Station, Tex. 

"Italicized number in parentheses refer to items in "Liter- 
ature Cited" preceding the appendix. 



1. Main progra'm. 

2. Subroutines. 

3. A data card containing the number of 
commands in the command table. 

4. A data, card containing the ZALFA array. 
B. Seventeen data cards containing the com- 
mand table. 

6. The users program deck consisting of 
program and data cards. 

A printout of the main program, subroutines, 
ZALFA array, and command table is given in 
the appendix. 
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The main program reads the command table 
and then calls the HONDO subroutine (2) to 
read a program data card. Subroutine HONDO 
determines the command name and number by 
comparing columns 1 through 20 of the program 
data card with the command table . Then 
HONDO determines individual data items by 
comparing columns 21 through 80 of the pro- 



gram data card with the ZALFA array. The data 
are placed in an array and returned to the main 
program. Based on the command number, the 
main program calls the proper subroutine to do 
the desired calculations. When the calculations 
are complete, control is returned to the main 
program, and HONDO is called again to read the 
next program card. 



The procedures used in HYMO were selected 
because of their accuracy, general applicability, 
practicality of inputs, and computational effi- 
ciency. For most watersheds the input is easy to 
obtain, and the procedures produce reasonably 
accurate results without excessive computer 
time. 

Hydrograph Computation 

When flood routing is performed, a watershed 
is divided into many small areas according to its 
hydraulic characteristics. The hydrographs from 
these areas must be estimated, since streamflow 
measurements are seldom available. A procedure 
for computing unit hydrographs was developed 
previously (4). A modification of this procedure 
is used in HYMO. Unit hydrographs are divided 
into three parts for computation (fig. 1). From 
the beginning of rise to the inflection point, ^o > 
the hydrograph is computed by the two- 
parameter gamma distribution equation 



(1) 




3456 
t/tp 



where q = flow rate in cubic feet per second jit time /, 

({u = peak flow rate in cubic feet per second, 

tp = time to peak in hours, 

and fi = dimensionless parameter. 

From tQ to ti (^ '- t Q + 2/0 the hydrograph 
is computed by the recession depletion equation 



K 



(2) 



where q$ = ftw l 'ate at the inflection point, 

^o = time at the inflection point, 
and K = recession constant in hours. 

From t\ to oo the recession depletion equa- 
tion becomes 



q = 



where 4j = flow rate at /!, 

and K\ = 3K = second recession constant. 



(3) 



The dimensionless shape parameter, n, is a 
function of K/t p , as shown in figure 2. The peak 
flow rate is computed by the equation 



BAQ 



(4) 



Figure 1. Dimensionless unit hydrograph. 



where B = a watershed parameter, a function of ;; as 

shown in figure 3, 

A - watershed area in square miles, 
and Q = volume of runoff in inches. 

Therefore, the entire unit hydrograph can be 

computed if K and t p are known. K and t p 

can -be determined by hydrograph analysis (4) 

for gaged watersheds. To compute K and t p for 

.ungaged watersheds, HYMO uses the equations 



K = 



and tp - 4.63/10.4225/^-0.46^/^0.133 (6 ) 

where SLP = difference in elevation in feet, divided by 
flood-plain distance in miles, between 
watershed outlet and most distant point 
on the watershed, 

and L/W - watershed length-width ratio. 

Storm hydrographs are computed by con- 
volving unit hydrographs with incremental 
source runoff. To compute incremental source 
runoff, the mass rainfall curve is broken into 
equal time increments, and the Soil Conserva- 
tion Service (SCS) rainfall -runoff relationship 
(J) is applied. The SCS rainfall-runoff relation- 
ship is expressed in a set of numbered curves. 
The SCS National Engineering Handbook (3) 
provides detailed instructions for selecting the 
proper curve number. 

Hydrographs computed by this procedure 
compared closely with measured hydrographs 
from 34 watersheds located in Texas, Oklahoma, 
Arkansas, Louisiana, Mississippi, and Tennessee. 
The watershed areas ranged from 0.5 to 25 
square miles. 

Flood Routing 
Streams and valleys 

The variable storage coefficient (VSC) flood- 
routing method (5) was selected for HYMO. 




Figure 2. Relationship between dimcnsionless shape 
parameter and recession constant/time to peak. 
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Figure 3. Relationship between dimensionless shape 
parameter n and watershed parameter B. 



The VSC method has been revised (6) to 
account for the variation in water surface slope 
during a flood. The revised VSC method is about 
as accurate as the implicit method (/) and has 
the general applicability of simpler storage 
methods. Although an iterative solution is used, 
the VSC method requires little computer time 
and is free of convergence problems. 
The VSC routing equations are 



0, = 
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In these equations subscripts 1 and 2 refer to 
the beginning and end of the time interval <v; the 
units are cubic feet per second for flow, hours 
for time, feet per second for velocity, and feet 
.for length and depth. The symbols are defined as 
follows: 

/ = inflow rate. 
= outflow rate. 

/1+/2 
I a = -) ~ average inflow rate, 

C = storage coefficient, 

7" = travel time through the reach. 

L = reach length. 

V = velocity. 

normal slope. 
D = depth. 



Since TI and 6*2 are dependent upon 02, an 
iterative technique is required to solve the 
routing equations. In equation 7, I a and are 
known, and C\ can be computed from equation 
9. This leaves only 2 and CT. as unknowns. 0\ 
can be used as a first approximation of O*i- The 
normal depth and velocity for the approximate 
value of 0-2 are entered into equation 11 for 
computing r 2 . Then equation 8 is used to 
compute GI. The second approximation of 0% is 
then obtained from equation 7. This iterative 
process continues until the difference between 
successive 2 values is acceptable. HYMO is set 
to accept differences of 0.1 percent or less. 
Usually about four iterations are required. 
Reservoirs 

HYMO uses the storage-indication method (3) 
to route floods through reservoirs. This method 
has been widely used and accepted because it is 
practical and accurate. The SCS National Engi- 
neering Handbook gives detailed instructions for 
using the method. 



Rating Curves 

Rating curves must be available at enough 
locations along a valley to adequately describe 
the hydraulics of the stream and valley. Most of 
these rating curves must be computed because 
there are never enough measured rating curves. 

HYMO uses Manning's equation to compute 
the normal flow-rating curves that are used in 
the VSC flood-routing method. The normal 
flood-plain slope is determined for each valley 
section by plotting a profile of the flood plain. 
The normal channel slope is determined by 
plotting a profile of the flood plain with channel 
distances. 

Sediment Yields 

The universal soil loss equation (8) was 
modified to compute the sediment yield for 
individual storms on watersheds. The modified 
equation is 

- 56 > 



X0 



(12) 



where 5 



K 

C 

P 

LS 



and 



sediment yield in tons, 

peak flow rate in cubic feet per second, 

volume of runoff in acre-feet, 

the soil-erodibility factor, 

the cropping-management factor, 

the erosion control practice factor, 

the slope length and gradient factor. 



Detailed instructions for determining K, C, P, 
and LS are given by Wischmeier and Smith (8). 

Since equation 12 was designed to compute 
sediment yield from watersheds, a delivery ratio 
is not needed. The delivery ratio is built into 
equation 12 by including the peak flow rate. 
Many of the watershed characteristics that in- 
fluence the peak flow rate also affect the 
delivery ratio. Equation 12 has performed well 
under limited testing, but future refinements are 
expected. 



RULES FOR USING HYMO 



The reader should refer to the table," "Ex- 
ample Input for HYMO Commands," as he 
follows the narrative description of the rules. 
The example HYMO program that is presented 
near the end of the manual may also be helpful. 

General Rules 

HYMO commands are expressed in the first 
20 columns of the data card, and columns 21 



.through 79 are used for numeric data and 
(keywords. Column 80 is reserved for a page, 
change code (an asterisk in column 80 causes the 
card to be printed on a new page). Continuation 
cards are allowed when 59 characters are insuf- 
ficient to express the data. 

The data can be written in any format, but at 
least one blank space must be left between data 
items. A decimal is required for numbers con- 



taining fractions, but not for whole numbers, may be used at any point in a HYMO program 
Keywords can be written with the data to by punching an asterisk in column 1 and the 
describe individual data items. Comment cards comment in columns 2 through 79. 



Example input for HYMO commands 



Command 



Required input 



START 
STORE HYD 

COMPUTE HYD 

PRINT HYD 

PUNCH HYD 

PLOT HYD 

ADD HYD 

STORE RATING CURVE 



RAINFALL BEGINS AT 12.5 HRS 



PUNCH CODE=1 



ID=1 HYDNO=301 DT-.2 HR DA=1.5SQMI 

FLOW RATES= 10 50 100 500 1000 1800 2000 1900 1500 1200 

1000 800 600 500 400 300 200 100 50 10 1 

ID=2 HYDNO=302 DT=.5 HR DA=2.1SQMI CN=90 

HT=100FT L=3.3MI 

MASS RAINFALL = .31 .61 1.04 1.84 2.74 3.06 3.45 4.33 4.75 

ID=2 CODE=1 



ID=2 



ID 1=3 ID 11=4 



ID=4 HYDNO<=101 ID 1=6 ID 11=6 



ID-2 VS NO=15 

ELEV AREA 

496.6 

497 

498 

499 



500 




2 
9 

19 
30 



FLOW 


1 

19 
52 
98 



COMPUTE RATING CURVE 



COMPUTE TRAVEL TIME 

ROUTE 

ROUTE RESERVOIR 



ERROR ANALYSIS 
SEDIMENT YIELD 
FINISH 



ID=1 VS NO=10 NO SEGS=3 MIN ELEV=482 FT MAX ELEV=492 FT 
CH SLP=.006 FP SLP-.0075 N=.OB DIST=175 FT 
N=-.03" DIST=205 FT N=.06 DIST=450 FT 

DIST ELEV DIST ELEV DIST ELEV DIST ELEV 
492.0 100 490.0 175 484.0 188 482.0 

190 482.0 205 484.0 250 486.0 275 488.0 

310 490.0 450 492.0 

ID=3 REACH NO=8 NO VS=5 L=4500 FT SLP=.0076 
ID=3 HYDNO=8 INFLOW ID=6 DT=.25 HR 

ID=5 HYDNO=501 INFLOW ID=1 

OUTFLOW (CFS) STORAGE (AC FT) 



22 533 

200 566 

1000 601 

2000 648 

3000 694 

ID 1=3 ID 11=6 

ID=6 SOIL-.34 CROP=.5 EP=.6 LS=.3 



Six hydrographs can be stored in a HYMO 
program at a time. The hydrographs are identi- 
fied by storage location numbers 1 through 6. 
Therefore, the same storage location number 
must be used for many hydrographs in a HYMO 
program. This is especially true when routing is 
done through large watersheds. However, no 
more than six hydrographs are ever needed at 
one time because HYMO programs begin at the 
head of a watershed and work downstream 
through one reach at a time. When a storage 
location number is used to store or compute 
another hydrograph, the first hydrograph is lost. 
The user should be sure that the hydrograph will 
not be referred to again before using the storage 
location number for another command. 

To store, compute, or route a hydrograph, the 
user must specify the time increment. There are 
no rigid rules about selecting the time incre- 
ment, but generally it should not be greater than 
one-fifth of the time to the peak of the 
hydrograph. This rule usually provides enough 
points to adequately define the hydrograph. All 
hydrographs are limited to 300 points. 

For the commands "STORE HYD," "COM- 
PUTE HYD," "ADD HYD," "ROUTE," and 
"ROUTE RESERVOIR," the user must specify 
the number of the outflow hydrograph. The 
hydrograph identification numbers are used to 
designate specific routing reaches, incremental 
areas, reservoirs, and partial hydrographs. The 
partial hydrograph number is given to all hydro- 
graphs other than outflow hydrographs from 
reaches, incremental areas, or reservoirs. The 
identification numbers for each group are 

Reaches 1-100 

Partial hydrographs 101-300 

Incremental areas 301-500 

Reservoirs 501+ 



Command Rules 

The first command for any watershed is 
START. The two data items associated with this 
command are the time rainfall begins on the 
watershed and a code for punching output data. 
If a storm is to be routed through a watershed 
only once, the punch code is deleted. However, 
if more than one routing is to be performed, set 
the punch code equal to a positive number, and 
the output data for the first routing will be 



punched for use in the second routing. More 
than one routing is usually required. 

Two commands, RECALL HYD and STORE 
TRAVEL TIME, were designed to be computer 
punched for second routings; consequently, 
these commands do not appear in the table. 

The STORE HYD command is used to store 
the coordinates of a hydrograph in the com- 
puter. It can be used for storing measured 
hydrographs or hydrographs computed by 
methods other than the one used in HYMO. The 
input data required for STORE HYD are storage 
location number, hydrograph identification 
number, time increment, watershed area, and 
flow rates of the hydrograph at the specified 
time increment. 

The COMPUTE HYD command is used to 
compute hydrographs from the incremental 
areas of the watershed. The first five items of 
data are storage location number, hydrograph 
identification number, time increment, water- 
shed area, and SCS runoff curve number (3), 
Normally, data items 6 and 7 are watershed 
height and main stem length. The height and 
length are used to compute the recession con- 
stant K and the time to peak f, } , However, if K 
and tp are known or estimated by some other 
method, they can be entered directly into the 
program. This is accomplished by placing a 
minus sign before the values of K. and t }) and 
entering them as data items 6 and 7, respective- 
ly. The remaining data items are values of the 
mass rainfall at the specified time increment. 

Since most watersheds have a limited number 
of rain gages, the same mass rainfall data may be 
used to develop several hydrographs. Once the 
mass rainfall data have been entered in a 
COMPUTE HYD command, they can be re- 
peated for any number of COMPUTE HYD 
commands without repunching the data. 
Instead, punch a negative number for the eighth 
data item of all COMPUTE HYD commands that 
use the same rain gage. When data from another 
rain gage are entered, the data from the first rain 
gage are lost and cannot be recalled by using the 
negative number code. 

The RECALL HYD command is one of the 
two commands that are computer punched. 
When the punch code is a positive number, the 
output from STORE HYD and COMPUTE HYD 
are punched on cards with the RECALL HYD 
command. The RECALL HYD command stores 



the computed and stored hydrographs on cards; 
it is therefore not necessary to recompute 
hydrographs for future routings. Instead, the 
previously computed hydrographs are read into 
the program, thus saving considerable computer 
time. 

Although the input data for the RECALL 
HYD command are never punched manually, a 
list of the data items may be helpful in checking 
computer-punched cards. The input data are 
storage location number, hydrograph identifica- 
tion number, time increment, drainage area, 
peak flow rate, runoff volume, number of 
hydrograph points, and flow rates of the hydro- 
graph. 

The PRINT HYD command is used to print 
coordinates of a hydrograph, volume of runoff, 
and peak flow rate . The required input data are 
the storage location number and a peak-volume 
code. The peak-volume code is deleted if a 
complete hydrograph printout is desired. If a 
printout of only the runoff volume and the peak 
flow rate is needed, the peak-volume code is set 
to a positive value. 

The PUNCH HYD command is used to punch 
any hydrograph in a HYMO program in the 
proper form for the RECALL HYD command, 
PUNCH HYD has two purposes: (1) If the 
punch code is not used, PUNCH HYD can be 
used to punch one or more hydrographs for 
future use; and (2) if it is desirable to punch 
outflow hydrographs associated with ROUTE, 
ROUTE RESERVOIR, or ADD HYD, PUNCH 
HYD must be used because the punch code only 
provides for punching hydrographs associated 
with STORE HYD and COMPUTE HYD. The 
only datum required for PUNCH HYD is the 
storage location number of the hydrograph to be 
punched . 

The PLOT HYD command is used to plot 
hydrographs in a HYMO program. It will plot 
one hydrograph on a set of axes, or if a 
comparison is desired, it will plot two hydro- 
graphs on the same set of axes. The required 
input data are the storage location numbers of 
the hydrographs to be plotted. 

The ADD HYD command adds the coordi- 
nates of any two hydrographs. The hydrographs 
are added at a time increment equal to that of 
the hydrograph with the shorter time increment. 
The only data required are the storage location 
number and hydrograph identification number 



of the added hydrograph and the storage loca- 
tion numbers of the two hydrographs to be 
added. 

The STORE RATING CURVE command is 
used to store rating curves that have been 
measured or computed previously. STORE 
RATING CURVE will save considerable com- 
puter time if measured or computed rating 
curves are available. The input data are the 
storage location number, valley, section number, 
and individual rating curve points described by 
elevation, end-area, and flow rate. The number 
of points used to describe a rating curve is 
limited to 20. 

The COMPUTE RATING CURVE command 
is used to compute the stage-area-flow relation- 
ship for a valley section. The input data are 
storage location number, valley section number, 
number of segments in the valley section, 
minimum elevation, maximum elevation, chan- 
nel and flood-plain slopes, Manning's n value and 
segment boundary point for each segment, and 
horizontal and vertical position of points de- 
scribing the valley section. 

The storage location numbers of the valley 
sections in a particular reach must begin with 1 
and increase by one for each valley section in 
the reach. However, the numbers are assigned 
without regard to upstream or downstream 
order. The valley section identification number 
can be any number from 0.1 to 999.9. These 
rules concerning storage location and valley 
section identification numbers also apply to the 
STORE RATING CURVE command. 

Normally, valley sections are divided into 
three segments (two flood-plain segments and a 
channel segment) for computing the rating 
curve. However, some valley sections may have 
more than one channel or may have an extreme 
variation in n values across the flood plain, thus 
requiring more than three segments, A maxi- 
mum of six segments is permitted. Manning's n 
values for each segment are input with seg- 
ment boundary point (distance from the begin- 
ning of the valley section to the end of the 
segment). Flood-plain n values are positive and 
channel n values are negative. 

Twenty points are used to define a rating 
curve. The location of the points is determined 
by dividing the difference between the maxi- 
mum and minimum elevations into 19 equal 
increments, 
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The COMPUTE TRAVEL TIME command is 
used to compute the normal flow travel time 
relationship used in ROUTE. The input data are 
storage location number, reach identification 
number, number of valley sections in the reach, 
reach length, and slope. The reach identification 
number can be any number from 0.1 to 999.9. 
The maximum number of valley sections per 
reach is six. The slope can be either the channel 
or flood-plain slope or a weighted average of the 
two. If flow is confined to the channel, the 
channel slope is of course applicable. If most of 
the flow is in the flood plain, usually the 
flood-plain slope is used. However, a weighted 
slope based on the relative rates of flow in the 
channel and the flood plain may be used. 

The COMPUTE TRAVEL TIME command 
considers each rating curve in the reach in 
computing the travel time flow relationship. 
COMPUTE TRAVEL TIME automatically 
selects the flow rates that are used in computing 
individual travel times. The flow rates of the 
rating curve with the lowest maximum flow rate 
are chosen. If the flow rates of any other rating 
curve in the reach were chosen, the rating curve 
with the lowest, maximum flow rate would have 
to be extrapolated. The travel time table is 
limited to 19 points because of the 20-point 
limit for rating curves. 

The STORE TRAVEL TIME command is one. 
of the two computer-punched commands. When 
the punch code is a positive number, the output 
from COMPUTE TRAVEL TIME is punched on 
cards with the STORE TRAVEL TIME com- 
.mand, Therefore, it is not necessary to re- 
compute rating curves or travel time for future 
routings. Instead, STORE TRAVEL TIME reads 
the previously computed travel time flow rela- 
tionship into the program, thus saving consider- 
able computer time. 

The input data for STORE TRAVEL TIME 
are not punched manually, but a list of data 
items may be helpful in checking computer- 
punched cards. The input data are storage 
location number, reach identification number, 



reach length, slope, and individual points of the 
relationship defined by depth, flow, and travel 
time . 

The ROUTE command is used to route floods 
through streams and valleys. The input data are 
storage location number and hydrograph identi- 
fication number of the outflow hydrograph, 
storage location number of the inflow hydro- 
graph, and time increment. The storage location 
number of the outflow hydrograph must be the 
same as the storage location number used in 
COMPUTE TRAVEL TIME for the reach. To 
prevent unnecessary program stoppage, ROUTE 
extrapolates the travel-time table when it is 
exceeded and writes the message, "TRAVEL 
TIME TABLE EXCEEDED." 

The ROUTE RESERVOIR command is used 
to route floods through reservoirs. The input 
data are storage location number and hydro- 
graph identification number of the outflow 
hydrograph , storage location number of the 
inflow hydrograph, and individual points of the 
reservoir's outflow-storage relationship. The out- 
flow-storage relationship must be expressed in 
20 points or less. If the outflow-storage relation- 
ship is exceeded, ROUTE RESERVOIR will 
extrapolate the relationship and write the mes- 
sage, "STORAGE-DISCHARGE TABLE EX- 
CEEDED." 

The ERROR ANALYSIS command is used to 
determine the error standard deviation and the 
percentage error in peak flow between any two 
hydrographs in a HYMO program. These func- 
tions make ERROR ANALYSIS useful in re- 
search. The input data are the storage location 
numbers of the two hydrographs to be analyzed, 
The SEDIMENT YIELD command is used to 
compute the sediment yield at any point in a 
watershed. Input data required are storage loca- 
tion number of the hydrograph from the area, a 
soils factor, a crop factor, a slope length and 
gradient factor, and a conservation practice 
factor (5). 

The FINISH command is used to end HYMO 
programs. There are no data associated with 
FINISH. 



EXAMPLE HYMO PROGRAM 



A short example problem is presented to 
demonstrate HYMO. Figure 4 is a map of the 
6,84-square-mile Brushy Creek watershed near 



Riesel, Tex. A flood will be routed ftirough'the 
watershed in its present condition, and the 
routed outflow hydrograph will be compared to 



the hydrograph measured at gaging station G. 
Also the sediment yield will be predicted and 
compared with the measured sediment yield. 
Then the same flood will be routed through the 
watershed with two proposed reservoirs. To 
determine the effects of the reservoirs, the 
outflow hydrograph and sediment yield will be 
compared to the outflow hydrograph and sedi- 
ment yield of the present-condition routing. 



Comment cards and keywords are used liber- 
ally in the example problem to acquaint the user 
with HYMO. After becoming familiar with 
HYMO, the user may write fewer comments and 
keywords, but generally users find them both 
quite helpful in describing the problem. To save 
space in the example problem, few of the 
hydrographs are printed or plotted. Some users 
may choose to print and plot all hydrographs. 



LEGEND 
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Figure 4, Brushy Creek watershed near Rlesel, Tex. 
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//$OPTIONS 

C HYHO A 1 

C A 2 

C THIS PROGRAM IS A PROBLEM ORIENTED COMPUTER LANGUAGE FOR A 3 

C DEVELOPING WATERSHED MODELS. A * 

C A 5 

C A 6 

C DEFINITION OF VARIABLES A 7 

C ID = STORAGE LOCATION NUMBER A 8 

C NHD - HYDROGRAPH IDENTIFICATION NUMBER A 9 

C DT = TIME INCREMENT IN HOURS A 10 

C DA = DRAINAGE AREA IN SQ MI A 11 

C OCFS = FLOW RATE ARRAY A 12 

C QCFSF = CUTOFF FLOW RATE IN CFS A 13 

C IENO = NUMBER OF POINTS IN A HYDROGRAPH A 14 

C SCFS - FLOW RATE ARRAY FOR TRAVEL TIME TABLE A 15 

C A = END AREA IN SQ FT A 16 

C Q = FLOW RATE ARRAY FOR RATING CURVES A IT 

C RAIN = MASS RAINFALL ARRAY A 18 

1 COMMON CFS(300),OCFS(300,6),IEND<6),DATA(310),DA<6),DP(20),NPU,NHD A 19 
1, SCFS120I,C<20) ,A{20,6),Q(20,6),RAIN(200) ,DEEP (20, 6 ) f NER, MAXNO ,NCO A 20 
2MM,ICC,NCODE,DIST(61,SEGN(6),CTBLE(50,11),ITBLE<50,2),ZALFA(20),DT A 21 
3(6),TIME,PEAK(6),ROIN,ISG(6) A 22 

2 NCODE=0 A 23 
C READ ZALFA ARRAY A 24- 

3 READ (5,20) ( ZALFAI I) , 1 = 1 , 1 5) A 25 
C READ NUMBER OF COMMANDS A 26 

4 READ 15,21) NCOMM A 27 
C READ COMMAND TABLE A 28 

5 READ (5,22) ( (CTBLE ( I , J ) , J=l , 11 J , ( I TBW I I , J ) , J= It 2 ) , 1=1,NCOMM t A 29 

6 WRITE (6,23) ( ZALFA ( I ) , 1=1 , 15) A 30 

7 WRITE (6,24) A 31 

8 WRITE 16,25) ( 1 CTBLE ( I , J ) , J = l , 11 ) , ( ITBLE ( I , J ) , J = l, 2) , 1 = 1 , NCOMM) A 32 
C ZERO CODES -- NPU * PUNCH, ICC CONTINUATION CARD, NER = ERROR A 33 

9 NPU=0 A 34 

10 ICC=0 A 35 

11 1 NER=0 A 36 

12 CALL HONDO A 37 

13 IF (NER) 2,2,19 A 38 

14 2 GO TO (3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, I7 t 18,19), NCODE A 39 

15 3 TIME=DATA(1) A 40 

16 NPU=OATA(2) A -VI 

17 GO TO 1 A 42 
IB 4 CALL STHYD A 43 

19 GO TO I A 44 

20 5 CALL RECHD A 45 

21 GO TO 1 A 46 

22 6 CALL CMPHYD A 47 

23 GO TO 1 A 4B 

24 7 CALL PRTHYD A 49 

25 GO TO 1 A 50 

26 8 CALL PUHYD A 51 

27 GO TO 1 A 52 

28 9 CALL HPLOT A 53 

29 GO TO 1 A 54 

30 10 CALL ADHYD A 55 

31 GO TO 1 A 56 

32 11 CALL SRC A 57 

33 GO TO 1 A 58 

34 12 CALL CMPRC A 59 
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35 




GO TO 1 


A 


60 


36 


13 


CALL SIT 


A 


61 


37 




GO TO 1 


A 


62 


38 


14 


CALL-CMPTT 


A 


63 


39 




GO TO I 


A 


64 


40 


15 


CALL ROUTE 


A 


65 


41 




GO TO 1, 


A 


66 


42 


16 


CALL RESVO 


A 


67 


43 




GO TO 1 


A 


68 


44 


17 


CALL ERROR 


A 


69 


45 




GO TO 1 


A 


70 


46 


18 


CALL SEDT 


A 


71 


47 




GO TO 1 


A 


72 


40 


19 


STOP 


A 


73 




C 




A 


74 


49 


20 


FORMAT 15AU 


A 


75 


50 


21 


FORMAT 12) 


A 


76 


51 


22 


FORMAT 2A1,9A2,2I3) 


A 


77 


52 


23 


FORMAT 1H1,9X,8HZALFA = ,15A1///) 


A 


78 


53 


24 


FORMAT 16X,13HCOMMAND TABLE//) 


A 


79 


54 


25 


FORMAT 10Xt2Al.9A2f2.I3) 


A 


80 


55 




END 


A 


81- 


56 




SUBROUTINE HONDO 


B 


1 




C 


THIS SUBROUTINE READS IN A DATA CARD , SEARCHES AN ALPHAMERIC 


B 


2 




C 


CODE TABLE TO DETERMINE THE NCODE OF THE OPERATION AND 


B 


3 




C 


COLLECTS VARIABLES FROM THE FREEFLOATING DATA FIELD 


B 


4 


5.7 




COMMON CF 5(300) ,QCFS< 300,6) , IEND (6 ) ,DATA1 310) , DA(6 ) , OP{ 20 I t NPU,NHD 


B 


5 






1,SCFS(20),C(20) ,A(20*6),Q(20,6) I RAIN(200),DEEP(20,6),NER,MAXNO|NCO 


B 


6 






2MM,ICC,NCODE,DIST(6),SEGN(6),CTBLE(50,11),IT8LE(50,2) ,ZALFA ( 20 ) , DT 


B 


7 






3(6),TIME,PEAK(6),RQIN,ISG(6) 


3 


8 


50 




DIMENSION CHAFU60), ALPHA(ll) 


B 


9 


59 




DIMENSION AUXA(10)i AUXB(IO) 


B 


10 


60 




IF ( ICC) lilt 3 


8 


11 




C 


READ IN DATA CARD 


B 


12 


61 


1 


READ (5,42) ( ALPHA{ I ) , 1 = 1 f 1 1 ) , t CHARt I) , I = U60) 


B 


13 




C 


IF FIRST CHARACTER IS BLANK THE CARD IS A CONTINUATION OF 


B 


14 




C 


PREVIOUS CARD. 


B 


15 


62 




IF (ALPHA( 1)-ZALFA< ID) 2,9,2 


B 


16 


63 


2 


IF (ICC) 3,3,40 


B 


17 




C 


ASTERISK IN COL. 80 MEANS SKIP TO NEW PAGE BEFORE PRINTING CARD 


B 


18 


64 


3 


IF (CHAR(60)-ZALFA(1D) 4,5,4 


B 


19 


65 


4 


WRITE (6,43) 


B 


20 


66 


5 


WRITE (6,44) (ALPHAU),I=l,ll),(CHAR(I)tI = lf60) 


B 


21 




C 


IF FIRST CHARACTER IS A * THE PREVIOUS CARD WAS A COMMENT CARD 


B 


22 


67 




IF (ALPHAC 1)-ZALFA(12J) 10,6,10 


B 


23 




C 


IF PUNCH CODE POSITIVE, COMMENT CARDS ARE PUNCHED. 


B 


24 


6B 


6 


IF (NPU) 8,8,7 


B 


25 


69 


7 


WRITE (7,45) (ALPHAU ), Il, 11), (CHAR(I) ,1 = 1,60) 


B 


26 


70 


8 


ICC=0 


B 


27 


71 




GO TO I 


B 


28 


72 


9 


WRITE (6,44) (ALPHA(l),I=l,lU,(CHAR(I),I = l,6Q) 


B 


29 


73 




GO TO 24 


B 


30 




C 


SEARCH FIRST TWO ALPHAMERIC CHARACTERS TO SEE IF THEY ARE NUMBERS 


B 


31 


74 


10 


ICC=l 


B 


32 


75 




DO 12 1=1,10 


B 


33 


76 




IF (ALPHA! ll-ZALFAtm 11,15,11 


6 


34 


77 


11 


IF (ALPHA(2)-ZALFA(I) 1 12,15,12 


B 


35 


78 


12 


CONTINUE 


B 


36 




C 


STATEMENT NUMBER 7 IS BRANCHED TQ IF NUMBERS ARE PRESENT 


B 


37 
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C [F NOT NUMBER SEARCH COMMAND TABLE FOR MATCH B 38 

C CALL FIRST 10 VALUES FROM PERMANENT DATA STORAGE B 39 

79 DO 14 1=1,50 I 7, 

80 DO 13 J=l,li " I' 

81 IF (CTBLEt I,J)-ALPHA{JM 14,13,14 ?J 
C SN 10=PART MATCH J 7^ 

82 13 CONTINUE a 44 
C IF THIS LOOP IS COMPLETED HE HAVE COMPLETE MATCH- CALL NCOOE B 45 
C AND MAX NUMBER AND EXIT LOOP [j J* 

83 NCODE=ITBLEU,U j' 

84 MAXNO=ITBLE(I,2) 7j 

85 GO TO 21 B 3^ 

86 14 CONTINUE * *; 
C IF MAJOR LOOPS FINISHED WITHOUT A MATCH WRITE ERROR MESSAGE B 51 

C AND SET NER = 1 B jj* 

87 NER=L B " 

88 WRITE (6,46) B JJ 

89 RETURN B " 
C CONVERT DIGIT INPUT CODE FROM ALPHAMERIC TO INTEGER FORM B 56 

90 15 NCQDE=G IT (ALPHA t Ii2i 1.1+0.5 B 57 
C FIND MAX NUMBER OF DATA ITEMS FOR THIS NCODE B 58 

91 DO 17 1=1,50 B 59 

92 IF (ITBLE(I,1)-NCODE) 17,16,17 B 60 

93 16 MAXNO*ITBLEUf2) B 6L 

94 GO TO 21 B 62 

95 17 CONTINUE 8 63 
C SEARCH DATA ROUTINE B 6/ t 
C SEE IF ANY DATA FOR THIS CARD B 65 

96 DO 19 1=1,50 B 66 

97 IF (ITBLEI I,1)-NCDDE) 19,18,19 B 67 

98 18 MAXNO=ITBLEU,2) B 68 

99 GO TO 20 B 69 

100 19 CONTINUE B 70 

101 20 CONTINUE B 71 

102 21 IF (MAXNO) 23,22,23 B 72 

103 22 RETURN B 73 
C ZERO ARRAYS AND COUNTERS B 74 

104 23 DO 47 1=1,310 B 75 

105 47 DATA ( 11=0, B 76 

106 NDATA=1 B 77 

107 24 NCHAR=0 B 78 

108 25 DO 26 1=1,10 B 79 

109 AUXA(1)=0. B 80 

110 26 AUXBU)=Q. B 81 

111 IT1=1 B 82 

112 IT2=1 B 83 

113 SIGN-l. B 84 

114 LOGIT=0 B 85 

115 KDGIT^O B 86 
C CARRY OUT DIGIT BY DIGIT SEARCH AND ACCUMULATION B 87 

116 27 NCHAR=NCHAR+1 B 88 
C HAVE WE CONSIDERED ALL CHARACTERS - RETURN IF SO B 89 

117 IF (NCHAR-60) 28,32,1 B 90 

118 28 DO 29 1=1,15 B 91 

119 IF (CHAR(NCHARI-ZALFAUM 29,30|29 B 92 

120 29 CONTINUE B 93 

121 GO TO 32 B 94 

122 30 GO TO 133, 33, 33, 33, 33, 33 t 33, 33, 33,33, 32, 27,36, 32, 31, 27), I B 95 
C SN 39 HANDLES SIGN CONTROL ON 1130 VERSION B 96 

123 31 SIGN UO B 97 
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124 GO TD 27 B 93 
C CHARACTER IS BLANK OR COMMA - DOES IT FOLLOW A DIGIT B 99 

125 32 GO TO (27,48) , III B 10 Q 
C CHARACTER IS A DIGIT - HAS A DECIMAL BEEN ENCOUNTERED B 101 

126 33 GO TD (34,35), IT2 B 102 

127 34 LDGIT=LDGim B 103 

128 ITU2 B 10 4 

129 AUXA(LDGIT)=CHAR(NCHAR) B 105 

130 GO TO 27 e 106 

131 35 KDGIT=KDGIT+1 8 IQ? 

132 MJXB(KDGIT)=CHAR(NCHARI B 108 

133 GO TO 27 B 109 
C CHARACTER IS A DECIMAL - DOES IT FOLLOW A DIGIT B 110 

134 36 GO TO (37, 3B), IT1 B 111 

135 37 IT1=2 8 112 

136 LDGIT=1 B n3 

137 38 IT2=2 B 114 

138 GO TO 27 B 115 
C ROUTINE TO CONVERT ALPHABETIC ARRAY TO FLOATING POINT NUMBER B 116 

139 48 DATA t NDATA ) = GI T{ AUXA , 1 ,LDGIT, 1 . H-GIT ( AUXQ , 1, 10,0. ) 8 117 

140 DATA INDATA=DATA(NDATA)*SIGN B 118 
C IS ALL DATA FURNISHED YES-RETURN NO INCREASE N DATA KEEP ON B 119 

141 IF (NDATA-MAXNO) 41,39,39 B 120 

142 39 ICC=0 B 121 

143 40 RETURN B 122 

144 41 NDATA=NDATA+1 B 123 

145 GO TO 25 B 124 
C B 125 

146 42 FORMAT (2A1 ,9 A2 ,60A 1 ) B 126 

147 43 FORMAT (1H1) -B 127 

148 44 FORMAT { 5X ,2A 1 , 9A2, 60A1 ) B 128 

149 45 FORMAT ( 2A 1 ,9A2 ,60A 1 ) 8 129 

150 46 FORMAT ( 10X,20HCOMMAND NOT IN TABLE) B 130 

151 END B 131- 

152 FUNCTION GIT ( TCARD , J , JLAST, SHI FT) C 1 

153 DIMENSION TCARD(IO), A(IO) C 2 

154 DATA A(l)/lHl/,Af 2J/1H2/,A(3)/1H3/, A(4)/1H4/ ( A(5I/1H5/,A(6)/H6/ C 3 

155 DATA A(7)/1H7/,A(8)/1H8/,A(9)/1H9/,A(10)/1HO/ C 4 

156 GIT=0. C 5 

157 TEN=10. C 6 

158 SUM=0. C 7 

159 DO 3 JNOW=J, JLAST C a 

160 TTEST=TCARD( JNOW) C 9 
C CHECK FOR LAST ENTRY C 10 

161 IF (TTEST.EQ.O.) GO TO 4 C 11 
C FIND NUMBER AND COMPUTE VALUE C 12 

162 DO 2 NUMB=1,10 C 13 

163 IF (TTEST-A(NUMBI) 2,1,2 C 14 

164 1 ZTEST=NUHB C 15 

165 IF (ZTEST.EQ.IO.) ZTEST=0. C 16 

166 ' SUM=SUM*TEN+ZTEST C 17 

167 GO TO 3 C 18 

168 2 CONTINUE C 19 

169 3 CONTINUE C 20 

170 4 IF (SHIFTI 6,5,6 C 21 

171 5 FI=JNOW-1 C 22 

172 SUM=SUM*(0.1*#FI) C 23 

173 6 GIT=SUM C 24 

174 RETURN C 25 
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175 END C 26 ~ 

176 SUBROUTINE STHYD D 1 
C THIS SUBROUTINE STORES THE COORDINATES OF HYDRQGRAPHS. 2 

177 COMMON CFS(300),OCFSI300,6),IEND16),OATA(310),DA<6),DP(201,NPU,NHD D 3 
1,SCFS(20),C(20) ,A<20,6) ,Q(20,6) , RA IN ( 200 I , DEEP ( 20, 6 ) ,NER, MAXNO.NCO D 4 
2MM,ICC,NCQDE,DIST(6),SEGN(6),CTBLE 150,11 I, ITBLF.150,2) , ZALFA ( 20 ) , DT D 5 
3(6). TIME, PEAKI6) ,ROIN,ISG(6J D 6 

178 ID*DATA(U D 7 

179 NHD=DATA{2) D 8 

180 DTUD)=DATA[31 D 9 

181 DAUD)=DATA<4) D 1 

182 J=5 D 11 
C REMAINING DATA ARE FLOW RATES D 12 

183 OCFS(1,ID)*DATA( J) D 13 

184 PEAKUD) = 1. D 1* 

185 RO=DATA(J D 15 
106 DO 4 1=2,300 D 1 & 

187 J-JH. D 17 

188 QCFS(I,ID)=DATA(JI D 18 

189 RQ=RQ+OCFS(I, ID) D 19 
C IS FLOW RECEDING 20 

190 IF (OCFSH,IOJ-OCFS(I-1IDH 1,2,2 D 21 
C HAS FLOW RECEDED TO CUTOFF RATE D 22 

191 1 IF (OCFSd ,IDM 5,5,4 D 23 
C DETERMINE PEAK FLOW D 24 

192 2 IF (OCFSII ,ID)-PEAKUD) ) 4,4,3 D 25 

193 3 PEAKUO)=OCFS( I ,ID1 26 

194 4 CONTINUE D 27 

195 5 IENDUD> = I-1 D 28 

196 M=IEND(ID> D 29 

197 RQ IN=1 RO*DT< I DM /(DA< 101*645.333 1 D 30 
C PUNCH CODE D 31 

198 IF (NPU) 7,7,6 D 32 

199 6 WRITE (7,81 I D,NHD, DT ( I D) ,DA( I ) ,PE AK ( I Dl , ROIN , I ENDl I D) D 33 

200 WRITE (7,9) ( OCFS ( J , 1 1 , J = l H ) D 34 

201 7 RETURN D 35 
C 36 

202 8 FORMAT! 'RECALL HYD ' , T2 1 , ' I 0= ' ,U,T29,'HYD NQ=', I 3 , T42 , ' OT= ' ,F9. D 37 

16, HRS' ,T61, 'DA=',F8.3 f SQ MI ' /T21 , ' PEAK* ' , F7. 0, ' CFS ' ,T40 , R0=' , D 38 

2F6.3, 1 INCHES', T59, 'NO PTS= , I 3/T2 1 , FLOW RATES') D 39 

203 9 FORMAT (T21,7F8.0) 40 

204 END D 41- 

205 SUBROUTINE RECHD F 1 
C THIS SUBROUTINE RECALLS PREVIOUSLY COMPUTED AND PUNCHED E 2 
C HYDROGRAPHS E 3 

206 COMMON CFS ( 300) ,DCFS1 300, 6) , I ENO< 6) ,DATA( 3 10 ) DAI 6 ) , DP< 20 ) , NPU NHD E 4 
1,SCFS(20),C(20) ,A(20,6) ,Q ( 20 ,6) , RAI N( 200 ) , DEEP ( 20 , 6 ) , NF.R , MAXNO, NCO E 5 
2MM,ICC,NCODEtOIST(6),SEGN(6),CTDLE(50,ll), ITDL E( 50 , 2 ) , ZALFA ( 20 ) , DT E 6 
3(6),TIME,PEAK(6),ROIN,ISG(6) E 7 

207 ID = DATA(1) E 8 

208 NHO=DATA(2) E 9 

209 DT(ID)=DATA(3) E 10 

210 DA(ID)=DATA(4) E 11 

211 PEAK(ID)=DATA(5) E 12 

212 ROIN=DATA(6) E 13 

213 IEND(ID)=DATA(7) E 14 

214 M=IEND(ID) E 15 

215 J=8 E 16 



216 

217 

218 I 
219 
220 



221 



222 



223 

224 
225 
226 
227 

C 

228 

229 1 
230 
231 

232 2 
233 
234 
235 
236 
237 

238 3 
239 
240 4 

C 

241 
242 
243 
244 
245 
246 
247 

248 5 

249 6 
250 

251 

C 

252 7 
253 
254 
255 
256 
257 

258 8 
259 
260 

261 9 
262 
263 
264 
265 



REMAINING DATA ARE FLOW RATES 

00 1 1 = 1, M 

OCFS( I ,ID)=DATA( JJ 

J = J+1 

RETURN 
END 

SUBROUTINE CHPHYD 

THIS PROGRAM DEVELOPS A UNIT HYDROGRAPH, CONVERTS MASS RAINFALL 
TO POINT RUNOFF, AND COMPUTES STORM HYDROGRAPHS BY SUMMATION. 
COMMON CFSOOO) , OCFSI 300, 6) , I END( 6 J ,DATA( 3 10 ) , DA( 6 I ,DP( 20 ) , NPU ,NHD 
1,SCFS[ 201,0(20) ,A(20,6) , Q ( 20 ,6) ,RAIN I 200) , DEEP { 20,6) ,NER , MAXNO, NCO 

2MM,ICC,NCODE,DISTl6),SEGNI&) t CTBLEI50,ll) T IT&LE(50,2),ZALFAI20l,OT 
3( 6) f TIME, PEAK! 6 ,ROIN,ISG(6) 

ID=DATA( 1) 

NHD=DATA(2) 

DT(ID)=DATA13) 

DAUDl=DATA(4) 

CN=DATA(5) 

ARE K AND TP FURNISHED OR WILL THEY BE COMPUTED 

IF (DATA(6M It2,2 

XK=-DATA(6) 

TP=-DATA{7) 

GO TO 3 

HT=DATAf 6) 

XL=DATA(7J 

SLOPE=HT/XL 

XLOW=(XL**2.) /DA{ ID) 

XK = 27.0*(DAUO)*#.23H*(SLQPE**(~,777))*(XLDW**.124) 

TP=4.63*{DA.( ID)*#,422)*(SLOPE**(-.46) ) * [ XLDW**. 133 I 

PEAK! IDI-U 

DO 4 1=1,300 

QCFSI I , 1U)=0. 

COMPUTE N BY ITERATION. 

XN=5.0 

XKTP=XK/TP 

DO 6 1=1,50 

TINFl.+SQRT{l./(XN-U) I 

XN1=.05/(XKTP*(ALOG(TINF/(TINF + .05) + .05)+l. 

DIFF=ABS(XM-XN) 

IF (DIFF-.001 ) 7,7,5 

XN=XN1 

CONTINUE 

WRITE (6,29) 

GO TO 28 

DETERMINE Cl. 

D6LT=TINF/100. 

TC1-0. 

XN1P=XN-1. 

XNlM^l.-XN 

DO B 1=2,101 

TC1=TC1+DELT 

CFS( I)=(TC1**XN1P)*EXP(XN1M*(TC1-1.)} 

SUM=CFS( IOD/2. 

DO 9 1=2,100 

SUN = SUM+CFS(I ) 

C1=SUM*DKLT 

CFSII=CFS(101) 

TT INF=TINF*TP 

TRECl=TTINF+2.*XK 



17 

ta 

L9 
20 
21 
22- 

l 
2 
3 

4 

5 

6 

7 

9 

9 
10 
11 
12 
13 
14 
15 
L6 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 

3n 

31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 
47 
4B 
49 
50 
51 
52 
53 



59 



E6E=EXP((TTINF-TRECn/XKI F 54 

267 XK1=3.*XK F 55 

268 B*645.333/(CH-CF$II*JXKTP*(l.-EEE)+EfiE*(XKl/TP)M F 56 
C COMPUTE B, QP, AND CFSI. F 57 

269 QP=(B*DA(IDn/TP F 58 

270 CFSI=QP*CFS(101) F 59 

271 CFSR1=CFSI*EEE F 60 

272 WRITE (6,30) XN.QP F 61 
C DETERMINE INCREMENTAL RUNOFF. F 62 

273 R=1000./CN-10. c 63 



III J=8 F 65 

276 IF (DATA(J)I 13,10,10 F 6& 

277 10 RAIN(1) = DATA(J) F 67 
HI D0 11 1 = 2,300 F 68 
279 J=J+1 F 6g 
290 RAINU = DATA( J> F 70 

281 IF (RAIN< M-RAINU-U) 12,11,11 F 71 

282 11 CONTINUE F 72 
233 U NUMB = I-1 p ,3 
284 13 DO 15 1 = 1, NUMB F 74 
295 IF (RAINm-Bl) 33,33,14 F 75 
286 33 DATA U 1 =0 . F 7! 
267 Qi^o. F 4; 

288 GO TO 15 F 73 

289 14 Q2 = i(RAlN(I)-Bl)**2.)/(RAINtlH-.8*R F 79 

290 DATA ( I) = Q2-Q1 F fifv 

291 Q1=Q2 J "V 

292 15 CONTINUE F gg 
C COMPUTE UNIT HYDROGRAPH. F B3 

293 T2=0. 



CFS(1J=0. F 85 

295 DO 20 1 = 2,300 F a6 

296 T2=T2fDTUD) F 87 

297 IF (T2-TTINF) 16,16,17 F 88 

298 16 CFS(I)=QP*((T2/TP)**XN1P)*EXP(XN1M*(T2/TP-1.)1 F 89 

299 GO TO 20 F 90 

300 17 IF (T2-TREC1) 18,18,19 F 91 

301 18 CFSm*CFSI*EXPTTINF-T21/XK) F 92 

302 GO TO 20 F ,3 

303 19 CFSm=CFSRl*EXPUTRECl-T2)/XKl) F 94 

304 IF (CFSm-1.) 21,21,20 F 05 

305 20 CONTINUE F 96 

306 1=300 F 07 

307 21 ICND=I F gg 

C COMPUTE STORM. HYDRQGRAPH. F 99 

308 DO 24 J=2,NUM8 F 100 

309 N = J + ICND-2 F JQI 

310 IF (N-300) 23,23,22 F 102 

311 22 N=300 

312 23 KK=J 



314 DO 24 K=KK,N F ,o 6 

315 aCFS<K,lD)=QCFS(K,IDH-DATA(J)*CFS<I) F 107 

316 24 1=1+1 p 1nn 

i*\-r M-U i F 10B 

317 H-K-1 F 

318 R0=0. p 

319 DO 26 1=2, M F 

C COMPUTE RUNOFF VOLUME F 112 

320 RO = RO+OCFS(I,IDI 



60 



C DETERMINE PEAK F U4 

321 IF (OCFSU , IDJ-PEAKUD)) 26,26,25 F 115 

322 25 PEAK(ID)=QCFS< I ,ID> F i, 6 

323 26 CONTINUE F U7 

324 I=M F ue 

325 IENDI ID) = I F U9 

326 M =I F 120 

327 ROIM = (RO*DT( IO)1/(DA(ID)#645.333) F121 



359 
360 



PUNCH CODE 



336 END 



343 GO TO 7 



346 GO TO 7 



F 



122 



328 IF (NPU) 28,28,27 F 12J 

329 27 WRITE (7,31) 1 D ,NHD ,DT( I D 1 , DAI I D > ,P EAK I I Dl ,ROI N, I ENDt ID I F 124 

330 HRITE (7,32) ( OCFS ( I , 10 J , 1= L,M ) F 175 

331 28 RETURN F U6 

C F 127 

332 29 FORMAT) 'N DID NOT CONVERGE AFTER 50 ITERATIONS.') F 12B 

333 30 FORMATIT10, 'SHAPE CONSTANT, N = ' , F6.3/T I0 r 'UNIT PEAK = ',F10.1t'C F 129 

1 1" J ft 

334 31 FORMAT! ' RECALL H YD ' , T2 I , ' 1 D = ' f I 1 ,T29, HYD NQ=',I 3 ,T42, *DT= ' ,F9. F 130 

16, HRS',T61,'DA='tFfl.3,' SQ HI ' /T21 , ' PgAK= ' , F7. a, 'CFS 1 ,T40 , ' R0=' , F 131 

2F6.3,' INCHES' ,T59, 'NO PTS= ' , I 3/T21 , * FLOW RATES'! F 132 

335 32 FORMAT (T21,7F8.0) F m 



337 SUBROUTINE PRTHYO G 1 
C THIS SUBROUTINE PRI NTS- THE COORD INA TES DF A HYDROGRAPH. G 2 

338 COMMON CFS ( 300 ) ,OCFSt 300 , 6) , I END (6 ) , DATA( 3 10 ) f DM 6 ] , DPI 20 ) ,NPU,NHD G 3 
L, SCFS<20) t C<203 , A ( 20, 6 J ,Q (20 , 6 ) , RAIN (200 ) , DEEP ( 20 , 6 I ,NER, MAXNQ.NCO G 4 

2MM,ICC,NCODE,DIST(6),SEGN(6I ( CTBLE (50,111 ,ITBLEt 50 ,Z) r ZALFA(ZO,OT G 5 

3(6^,TIHF. 1 PEAK{6I ,ROIN,ISG(6I G 6 

339 IO=DATAI1) G 7 

340 NPK=DATA(2 ) G 8 
C DETERMINE TYPE OF HYDROGRAPH G g 

341 IF (NHD-100) 6,6,2 G , 

342 1 WRITE (6,14) NHD G u 



G 12 



344 2 IF (NHD-300) 3 T 3,4 G ,3. 

345 3 WRITE (6,15) NHD G 



G 15 



347 4 IF (NHD-500) 1,1,5 G 16 

34B 5 WRITE (6,16) NHD G 17 

349 GO TO 7 G , e 

350 6 WR ITE (6,17) NHD G 19 
C POSITIVE NPK MEANS PRINT ONLY PEAK AND VOLUME G 20 

351 7 IF (NPK) 8,8,11 G 21 

352 8 J=0 G 22 

353 WR ITE (6,18) G 23 

354 M=IEND( ID) G ?4 

355 TIMEUTIME G 2 ~5 
C BUILD TIME ARRAY IN DATA G 2 6 

356 DO 9 1*1, M G 2 7 

357 DATA { I) = TIME1 G 2Q 
356 9 TIME1=TIME1+DT( ID) G 29 



G 30 

G 31 



361 10 J-J+1 G 32 

362 WRITE (6,19 I DAT A( U , OCFS ( I , ID t , 1 = J ,H , H5) G 33 

363 IF (J-M5) 10,11,11 G 34 

364 11 WRITE (6,20) RO I N f PEAK ( I D ) G 35 

365 IF <NPU) 13,13,12 G 36 

366 12 WRITE (7,21 ID, NPK G 37 
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367 


13 


RETURN 


G 


38 




C 




G 


39 


368 


14 


FORMAT (1HO,46X,21HHYDRGGRAPH FROM AREA ,13/1 


G 


40 


369 


15 


FORMAT (1HO,41X,19HPARTIAL HYDROGRAPH ,I4/) 


G 


41 


370 


16 


FORMAT (1HO,39X,29HOUTFLOH HYDROGRAPH RESERVOIR ,I4/) 


G 


4Z 


371 


17 


FORMAT (lHO f 44X,25HQUTFLOW HYDROGRAPH REACH ,I4/ 


G 


43 


372 


18 


FORMAT (10X,114HTIME FLOW TIME FLOW 


G 


44 






1TIME FLOW TIME FLOW TIME FLOW/1 


G 


45 






21X,113HHRS CFS HRS CFS HRS 


G 


46 






3 CFS HRS CFS HRS CFS) 


G 


47 


373 


19 


FORMAT |5(5X,F10.3,F10.0) 1 


G 


48 


374 


20 


FORMAT ( 1HO,9X,16HRUNOFF VOLUME = F10.3.BH INCHES /10X,22HPEAK OlS 


G 


49 






ICHARGE RATE = ,F10.L,4H CFS///) 


G 


50 


375 


21 


FORMAT! 'PRINT HYD' , T21 ' ID^ 1 t II T29, * CODE= ' , I 1 ) 


G 


51 


376 




END 


G 


52- 


377 




SUBROUTINE PUHYD 


H 


l 




C 


THIS SUBROUTINE PUNCHES HYDROGRAPHS IN FORM TQ BE USED BY 


H 


2 




C 


SUBROUTINE RECHD 


H 


3 


378 




COMMON CFS (300) ,OCFS(300,6I , I END 16) , DATAl 3 10) , DA ( 6 ) , DP { 20 > ,NPU , NHD 


H 


4 






liSCFSI 20),C(20) ,A(20,6) ,0(20,6) , RAI N(200 > , DEEP { 20 , 6 ) , NER, MAXNO,NCD 


H 


5 






2MM,ICC,NCODE,DIST(6),SEGN(6) ,CTBLE (50 , 1 1 ) , I TBLE < 50, 2 , ZALFA ( 20 t , DT 


H 


6 






3(6),TIME,PEAK(6),RQIN,ISG(6) 


H 


7 


379 




ID-DATA(l) 


H 


8 


380 




M=IEND( ID) 


H 


9 


381 




WRITE (7,1) ID,NHD,DT( IDI,DA{ ID) tPEAKUDI ,ROIN t lENDUD) 


H 


10 


382 




WRITE (7,2) (OCFSUtIDItI = L,M) 


H 


11 


383 




RETURN 


H 


12 




C 




H 


13 


384 


1 


FORMATf 'RECALL HYD ' ,T21 , ' I 0=' , U ,T29, 'HYD NQ= ' , I 3 ,T42 , ' DT= ' ? F9. 


H 


14 






16,' HRS' ,T61 , 'DA=' ,F8.3,' SQ Ml ' /T2L , ' PEAK= ' , F7.0 , 'CFS ' ,T40 , ' RQ=' , 


H 


15 






2F6.3,' INCHESST59,'NQ PTS= , I3/T2 1 , ' FLOW RATES') 


H 


16 


385 


2 


FORMAT (T21,7Ffl.O) 


H 


17 


386 




END 


H 


18- 


387 




SUBROUTINE HPLOT 


I 


1 




C 


THIS SUBROUTINE PLOTS EITHER 1 OR 2 HYDROGRAPHS ON A SET OF AXIS 


I 


2 


388 




COMMON CFSi300),OCFS(300,6t , I ENDt 6 > , DATA ( 3 10 ) , DA( 6) , DP( 20 ) , NPU ,NHD 


I 


3 






1, SCFS(20),C(20) ,A(20,6) ,0(20,6) , RAIN (200) , DEEP (20, 6) ,NER, MAXNO,NCO 


I 


4 






2MM,ICC,NCOD,DIST(6t,SGN(6),CTBLE(5Q,ll),!TBLE(50,2) , ZALFA ( 20 ) ,DT 


I 


5 






3(6),TIME,PEAK(6),ROIN,ISG(6) 


I 


6 


389 




ID1=DATA(1) 


I 


7 


390 




ID2-DATAI2) 


I 


a 


391 




DATA ZERO, PLUS, BLANK, DASH, DOT/ '0 ' , + ' , ' ', '-','./ 


I 


9 


392 




MRTO=1 


I 


10 


393 




XMRTO=1. 


I 


11 


394 




MAX=118 


I 


12 


395 




J = l 


I 


13 




C 


ARE THERE 1 OR 2 HYDROGRAPHS 


! 


14 


396 




IF IID21 1,1,2 


I 


15 




C 


DETERMINE HIGHEST PEAK IF 2 HYDROGRAPHS 


I 


16 


397 


1 


QMAX=PEAK( ID1) 


1 


17 


398 




GO TQ 14 


I 


18 


399 


2 


IF (PEAK(ID1)-PEAK(ID2) } 3,3,4 


I 


19 


400 


3 


QMAX=PEAK( ID2) 


I 


20 


401 




GO TO 5 


1 


21 


402 


4 


QMAX=PEAK( 101) 


I 


22 




C 


IF 2 HYDROGRAPHS DETERMINE LARGEST DT AND INTERPOLATE OTHER 


I 


23 




C 


HYDROGRAPH IF NECESSARY 


I 


24 


403 


5 


IF (DTUOlJ-DTt 1021 I 6,13,7 


I 


25 
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404 6 L=I01 T 

405 K=ID2 37 

406 GO TO 8 ;' 

407 7 L=ID2' 1 ;J 

408 K=ID1 ** 

409 8 M=IEND(LI [ ,r 

410 TID=DT(K) JJ 

411 TIDH=0. } " 

412 00 il 1*2, M ! " 

413 TIDH = TIDhH-DT(U j fj 

414 IF (TID-TIOH) 10,9,11 ,? 

415 9 J = J + 1 [ " 

416 CFS(J)=OCFS(I ,L) ,' 

417 T ID=TID*OT(K) 

418 GO TO 11 " 

419 10 J=J+1 j ^i 

*2 CFS(JIOCFS(I-l,LH-(aiO-TIOH+OT(LII/OT(L)*|OCFStI,L)-nCFS(I-l,L) I 42 

421 TIC=TIO+DT(K) } ^ 

422 11 CONTINUE J Tt 

423 IEND(L)=J { 7^ 

424 DT(L)=DT(K j 7* 

425 DO 12 1=2, J { ' 

426 12 UCFS(I,Lt=CFS(I J 7 

427 13 IF ( IEND( ID1I-IEN01 1021 ) 14,14,15 , *n 

428 14 M=IENP(ID1) '.' 

429 GO TO 16 

430 15 M=IEND(ID2) 

431 16 IF (M-MAX) 17,17,18 
C DETERMINE TIME SCALE 

432 17 MRTO=MAX/M 

433 XMRTQ=MRTO 

434 GO TO 19 

435 18 M=MAX 

436 19 YSCL=QMAX/50. 

C PLOT HYDRQGRAPHS ., 

437 DO 20 1=1, MAX 

438 20 CFS(I)=DASH 

439 WRITE (6,41) QMAX, ( CFS ( U , I = 1 ,MAX) , DOT 

440 Ql^QMAX 

441 Jl=10 

442 DO 37 J=l,50 

443 IF (J-Jli 23,21,23 

444 21 DO 22 1=1, MAX 

445 22 CFS( I)=DASH 

446 GO TO 25 

447 23 DO 24 1 = 1, MAX 

448 24 CFS(I)=BLANK 

449 25 Q2=Q1-YSCL 

450 K=l 5 

451 DO 28 1=2, M 7 , 

452 K=K+MRTO f 

453 IF (OCFSt I ,IDl)-Ql) 26,27,26 r 78 

454 26 IF (OCFSU ,ID1)-Q2) 20,28,27 70 

455 27 CFS(K) = ZERO 'f. 

456 28 CONTINUE oV 
^57 WRITE (6,44) DDT, (CFS ( I) , 1= l.MAX ) ,DOT R? 
^58 IF (ID21 34,34,29 
^59 29 K-l " 
460 DO 33 1=2, M 
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461 K = K+MRTO J 86 

462 IF (OCFS(IiID2)-QU 30,31,32 J J' 

463 30 IF (OCFS( I,ID2)-Q2) 32,32,31 { 88 

464 31 CFS(K=PLUS J * 

465 GO TO 33 J 

466 32 CFS(K)=BLANK J 9 J 

467 33 CONTINUE J ^J 

468 WRITE (6,42) ICFSU) , I1,MAX) ! 93 

469 34 IF (J-J1) 36,35,36 * 



470 35 Jl=JltlO 



1 95 



471 WRITE (6,43) Q2 ! 96 

472 36 Q1=Q2 J " 

473 37 CONTINUE l 9fl 

474 CFS(U=TIME l " 

475 DTT=DT( ID 1 ) *10./XMRTQ l 1 
C PUT TIME ARRAY IN CFS AND WRITE TIME SCALE I 101 

476 DO 38 1=2,12 J 1 2 

477 38 CFSI U = CFS( I-n+DTT I 103 

478 WRITE (6,45) (CF S( I J, 1=1 , 12 ) I 10* 

479 WRITE (6,46) * 105 

480 IF (NPU) 40,40,39 * 106 

481 39 WRITE ( 7, 47 - IDl , ID2 * 107 

482 40 RETURN l 10a 
C I 109 

483 41 FQRMAT(1X,F6.0,' CFS.',119A1) I HO 

484 42 FORMAT ( IH+ , 1 IX , 11BA1 ) * 111 

485 43 FORMAT {1H + .F6.0I I HZ 

486 44 FORMAT (11X,120A1I I 113 

487 45 FORMAT (6X,12F10.3 I H* 

488 46 FDRMAT(49X,'TIME HOURS'///) I 115 

489 47 FORMAT! 'PLOT HYD',T21,MD 1 = ' , II ,T29, ' I D 11 = ', ID I 116 

490 END I 



491 SUBROUTINE ADHYD J I 
C THIS SUBROUTINE ADDS TWO HYDROGRAPHS. J 2 

492 COMMON CFS (300) ,OCFS( 300, 6} , TEND ( 6) ,DATA( 310) , DAI 6 ) , DP( 20 ) , NPU ,NHD J 3 
1,SCFS(20),C(20 ,A{20 t 6),Q(20,6),RAIN(200) , DEEP ( 20 , 6 ) , NER ,MAXNO ,NCO J 4 
2MM ( ICC,NCOOE,OIST16),SEGN(6),CTBLE150,H),ITBLE<50,2),ZALFA<20),OT J 5 
3(6),TIME,PEAK(6),ROIN t ISG(6) J 6 

493 ID=DATA11) J 7 

494 NHD=DATA(2) J 8 

495 ID1=DATAU) J 9 

496 ID2=DATA(4) J 10 

497 PEAKUO) = 1. J U 
C MAKE TIME INCREMENTS EQUAL IF NOT EQUAL. USE SMALLER INCREMENT J 12 

498 IF (DT(ID1)-DT(ID2) ) 1,3,2 J 13 

499 1 OTUO)=OT( IDl ) J 14 

500 L=ID1 J 15 

501 K=ID2 J 16 

502 GO TO 6 J 17 

503 2 DT(ID=OTU02) J 10 

504 L=ID2 J 19 

505 K=I01 J 20 

506 GO TO 6 J 21 

507 3 OT(ID)=DT( IDl) J 22 

508 IF UENDUD1)-IEND(ID2)I 4,4,5 J 23 

509 4 M3 = IENDUDU J 24 

510 KUID2 J 25 

511 IEND(ID) = IENDUD2) J 26 

512 GO TO 18 J 27 
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5 M3 = IENDUD2) 

514 K1=ID1 J 28 

515 IEND< ID)=IENOUD1) J 29 

516 GO TO 18 J 30 
C DETERMINE DURATIONS OF FLOW J 3l 

517 6 XIEND1=IENDU01)-1 J 32 

518 XIEN02=IENDUD2)-1 ^ 33 

519 DUR1=X IEND1*DT( ID1I J 34 

520 DUR2=XIEND2*DT( ID2) > 35 

521 IF (DUR1-DUR2) 7,8,8 J 36 

522 7 IENDI ID)=DUR2/DT( ID 1 + 1. J 37 

523 M3=DURl/OT(IOI+l.- J 38 

524 K1=ID2 J 39 

525 GO TO 9 J 40 

526 8 I END! ID)=OUR1/DT( ID 1 + 1. J 41 

527 M3=DUR2/DT(ID+1. -J 42 

528 Kl-IDl J 43 

529 9 IF ( IEND{ ID)-300I 11,11,10 J ^ 

530 10 IEND(ID)=300 J 45 

531 11 M2=IENO(K) J 46 

532 J=l J 47 

C INTERPOLATE ONE HYDRQGRAPH IF NECESSARY J 48 

533 TIDH=0. J 49 

534 TID = DT( ID) J 50 

535 00 15 1=2, M2 J 51 

536 TIOH=TIDH+DT(K) J 52 

537 12 IF (TIDH-TID) 15,13,14 J 53 

538 13 J=J+1 J 54 

539 DATA (J)=OCFS(I,KI J 55 

540 TID = TID + DT( ID) J 56 

541 IF IJ-300) 15,16,16 J 57 

542 14 J=J+1 J 58 

543 J3ATA (J ) =OCFS ( I - 1 ,K M ( T I D-T IDH+DT. K ./DT ( K 1 * (OCFSU , K.-QCFS ( 1-1 , J 60 

544 TID=TIO + DT( ID I J 6l 

545 IF (J-300) 12,16,16 J 62 

546 15 CONTINUE J 63 

547 16 IEND(K)=J J ^ 

548 DO 17 1=2, J J 65 

549 17 OCFSU, K)=DATA(I| J 66 

550 18 M=!END< ID) J 67 

551 DA(ID) = DA( ID1I+OA(ID2) J 6S 

552 R0=0. J 69 
C ADD HYDROGRAPHS J 70 

553 DO 20 1=1, H3 J 7l 

554 OCFS(I,ID=QCFS(I,ID1)+OCFS(I,ID2) J , l\ 

555 IF (OCFSU, IDI-PEAKUDM 20,20,19 , If 

556 19 PEflK(ID)=OCFS(I,IOI , ll 

557 20 RO=RO+OCFS1 I, ID) *J 75 

558 IF (PEAK(ID-PEAK(K1H 21,22,22 , 77 

559 21 PEAKt ID)=PEAK(K1) , 'I 

560 22 IF (M-M3J 25,25,23 J . 1 

561 23 M3=M3+1 *> 79 

562 DO 24 I=M3,M , ^ 

563 OCFS(I,ID)=OCFS(I,K1I , JJ 

564 24 RO=RQ+OCFS(I, ID) , J, 

565 25 RQIN=<RQ*DTUDM/(OA(ID>*645.333) , H 

566 IF (NPU) 27,27,26 J . JJ 

567 26 WRITE (7 r 28) I D ,NHD , 1 01 , 1 02 1 Jl 

568 27 RETURN ^ 

J 87 



65 



C J 88 

569 28 FORMAT! 'ADD HYD' ,T21 , ' I D= ', 1 1 ,T2 9, HYD NO* , 13,145, ' ID I=',I1, J 89 

1T60,'ID 11=' ,11) J 90 

570 END J 91- 

571 SUBROUTINE SRC K * 
C THIS SUBROUTINE STORES AN ELEVATION - END AREA - FLOW TABLE. K 2 

57Z COMMON CFS (300 ) ,OCFS< 300, 6) , I END t 6 ) , DATA< 3 10) , DA( 6 ) , DP( 20 f NPU, NHD K 3 

1,SCFS(20),CI20),A{20,6) ,Q (20,6 ) , RAIN < 200 ) f DEEP I 30 , 6 ) , NER, MAXNO,NCO K 4 

2MM,ICC,NCaDE,DIST<6),SEGN(6),CTBLE(50,ll), ITBL E ( 50 , 2 ) , ZALF A ( 20 ) , DT K 5 

3(6), TIME, PEAK(6},RQIN,ISG(6} K 6 

573 IO=OATAtl) K 7 

574 VS=DATA(2 * B 
C VALLEY SECTION NUMBER K 9 
C REMAINING DATA ARE ELEVATION, AREA, AND FLOW FOR EACH POINT OF K 10 
C THE RATING CURVE K 11 

575 EMIN=DATA(3) K 12 

576 J = 3 K 13 

577 DO 1 1=1,20 K 14 
57B ELEV=DATA( J) K 15 

579 DEEP! I , I ) =DATA t Jl-EM IN K 16 

580 A( I ,ID)=DATA( J-t-U K 17 

581 Q! I,ID=DATA< J+2 I K 18 

582 J=J+3 K 19 

583 I CONTINUE K 20 

584 RETURN K 21 

585 END K 22- 

586 SUBROUTINE CMPHC L 1 
C THIS SUBROUTINE COMPUTES THE DISCHARGE END-AREA ELEVATION L 2 

c RELATIONSHIP FOR A VALLEY SECTION. i 3 

587 COMMON CFS ( 300) ,OCFS ( 300, 6) , I END (6 ) ,DATA( 3 10) , DA ( 6 I , DP ( 20 ) , NPU, NHD L 4 
l.SCFSI 20) ,CI20) ,A(20,6) ,0(20,6) ,RAIN(200) , DEEP (20, 6) ,NER , MAXNO, NCO L 5 
2MM,ICC,NCODE,DIST(6),SEGN{6),CTBLE(50,11), ITBL E ( 50 , 2 ) , ZALFA ( 20 , DT L 6 
3(6), TIME, PEAK(6>,RQIN,ISG<6) L 7 

588 1D=DATA(1) L 8 
C STORAGE LOCATION NUMBER. (1-6) L 9 

589 VS=OATAI2> L 10 
C VALLEY SECTION IDENTIFICATION NUMBER. L 11 

590 NSEG=DATA<3) L L2 
C NUMBER OF SEGMENTS IN THE VALLEY SECTION. L 13 

591 ELODATAU) L 14 

592 EMAX=DATA(5) L 15 
C MAXIMUM ELEVATION FOR COMPUTATIONS. L 16 

593 SLOPEl=DATA(6) L 17 
C CHANNEL SLOPE. L IS 

594 SLOPE2=DATA(7) L 19 
C FLOOD PLAIN SLOPE. L 20 

595 DIF=(EMAX-ELOt/19. L 21 

596 C( 1)=ELO L 2 2 

597 DO 1 1=2,20 L 23 

598 1 C{ I)=C( 1-1 J+DIF L ? ^ 
C SET AREA AND DISCHARGE ARRAYS =0. L 25 

599 DO 2 1=1,20 L 2 6 

600 A( 1,101=0. I ?7 

601 2 Q( I,ID)=0. L 28 



602 J=8 



L 29 



603 WRITE (6,24) VS L 30 
C READ N VALUES AND SEGMENT BORDER POINTS. L 31 

604 DO 3 I=1,NSEG L 32 



66 



605 SEGN1 I )=DATA( J J 

606 DISTI I)=DATA( J+l) , 

607 3 J = J + 2 L 



608 

609 DO 6 IMiNSEG i 37 

610 4 J = J + 2 L 3g 

AH IF (DATAl J)-OIST(II) 4,5,5 L 39 

612 5 ISG( I)=J+l L 4 

613 6 CONTINUE L 41 

C COMPUTE DISCHARGES AND END AREAS FDR EACH SFFMPMT L ^ 

6t4 DO 22 K=i,NSEG SEGMENT. L ^ 3 

615 J=JJJ L 4.4 

616 JJJl=JJJfl L 45 

617 IF (SEGN(KJ) 7,7,8 L 46 

618 7 SLOPE=S1_OPE1 L 47 

619 SEGN(K)=-SEGN(K) L 48 

620 GO TO 9 L 49 

621 8 $LCPE=SLDPE2 L 50 

622 9 SLPN=l.486*SLOPE**. 5 L 51 

C COMPUTE AREA AND DISCHARGE FOR SEGMENT L 52 

623 DO 21 1-2,20 scbMbHl. L 53 

624 AA=0. i 54 

625 P = 0. l_ 55 

626 J=JJJ-1 L 56 

627 DEP2=0. L 57 

628 10 J=J+2 L 58 

629 IF (J-ISG(KI) 12,12,11 L 59 
6 3 11 IP (AA-.001) 21,21,20 L 60 

631 12 IF (DATA(J)-CU I I 13,10,10 L 61 

632 13 OEP1=C( M-DATAI J) L 62 

633 IF (J-JJJ1) 16 t 16,14 L 63 

634 14 XL=DATA( J- 1 I-DAT A ( J-3 ) L 64 

6 35 DEP3 = ABS(OATA(J-2)-DATA<JM L 65 

636 XL=XL*DEP1/DEP3 L 66 

637 15 AA = AA+XL*( OEPH-DEP2J/2. L 67 

<> 38 P=P+SQRT((DEP1-DI:P2)**2*XL**21 L 68 

639 16 DEP2*DEP1 L 69 

640 J=J+2 L 70 

6 41 IF (J-ISG(K)) 17,17,20 L 71 
6 *2 17 IF (DATA(J)-C( [ H 18,10.19 L ^2 

643 18 DEP1=C{ I J-DATAt J) L 73 

644 XL=OATA< J- 1 J-OATA { J-3 ) L 74 

645 GO TO 15 L 75 
M> 19 OEPUO. L 76 

647 XL=()ATA(J-l)-DATA<J-3) L 77 

648 DEP3=ARS(DATA( J-2)-DATA(JM L 
6*9 XL=XL*DEP2/DEP3 L 79 

650 AA = AA + XL*(DEPl-t-OtEP2)/2. L 80 

651 P = P+SQRT< (DEP1-DEP2*'*2+XL**PI L 1 

652 DEP2=0. L 62 

653 GO TO 10 L R3 

654 20 K = AA/P L B4 

655 SGN=SEGN(K J-.0025*R L 85 

c 5aEv S s??aN AN " AREAS FOR ALL SEGMENTS ""IN TOTALS FOR t S? 



L 

658 21 CONTINUE L 90 

659 JJJ=J-3 L 91 

L 92 
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660 22 CONTINUE 

661 DO 23 1 = 1,20 L 

662 DEEPU,IO)=C(I)-ELO L 

663 WRITE (6,25) C t II , A ( I , ID) ,Ql 1 t I D ) L 



664 23 CONTINUE 



L 97 



665 RETURN L 98 
C L 99 

666 24 FORMAT(1HO,T42, 'RATING CURVE VALLEY SECT ION* , F5 . 1/T46 ,' WATER ', T56, L 100 

1'FLOW' , T66, 'F LOW' /T45 , ' SURF ACE 1 ,T56,'AREA' ,T66 , 'RATE ' /T46 , ' ELE V ,T L L01 

256, 'SQ FT' ,T66,'CFS' I L l02 

667 25 FORMAT (40X ,F 10. 2 ,2F10. 1) L 103 

668 END L L0/t - 

669 SUBROUTINE STT H l 
C THIS SUBROUTINE STORES A DEPTH - FLOW - TRAVEL TIME TABLE. M 2 

670 COMMON CFS ( 300} ,OCFS( 300,6) , I END ( 6 ) ,DATA< 3 10) , DAI 6 ) , DP< 20 ) , NPU ,NHD M 3 
1,SCFS[20),C(20) , At 20, 6) ,Q (20 ,6 } , RAINl 200 > , DEEP ( 20 , 6 ) ,NER, MAXNO, NCO M 4 
2MM,ICC,NCODE,DIST(6),SEGN(6},CTBLE<50,11) , ITBL Et 50 , 2 I , ZALF A( 20 ) , DT M 5 
3<6),TIME,PEAK(6),ROIN,ISG(6) H 6 

671 ID=DATA(1) H 7 

672 REACH=DATA(2) M 8 

673 XL=DATA(3) M 9 

674 SLQPE=DATA(4 H 10 

675 DIST(ID)=SLOPE*XL M 11 

676 J=5 M 12 

677 00 1 1=1,19 M 13 

678 OP(I) = DATA( J) M 14 

679 SCFSU )=DATA( J+1J M 15 

680 C( II-DATAt Jt2l H 16 

681 1 J=J+3 M 17 

682 RETURN M 18 

683 END M 19- 

684 SUBROUTINE CMPTT N 1 
C THIS SUBROUTINE COMPUTES THE TRAVEL TIME AT GIVEN N 2 
C DISCHARGE RATES N 3 

685 COMMON CFS ( 300) ,OCFS( 300,6) , IEND{6) ,DATA{3 10 ) , DA( 6 ) ,DP ( 20 ) ,NPU ,NHD N 4 
1,SCFS(20),C(20) ,AI20,6),Q(20,6) , RAIN ( 200 ), DEEP ( 20 , 6 ), NER, MAXNO, NCO N 5 
2MM,ICC,NCQDE,DIST(6),SEGN(6) ,CTBLE ( 50 , 1 1 ) ITBLE(50,2) ,ZALFA( 20 ) , DT N 6 
3{6),TIME,PEAK(6),ROIN,ISG<6) N 7 

686 ID=DATA(1) N S 
6B7 REACH=DATA(2) N 9 

688 NQVS=DATA(3) N 10 

689 XL=DATA14) N 11 

690 SLOPE=DATA<5) N 12 

691 DIST(ID)=SLOPE*XL N 13 

692 XLD36=XL/3600. N 14 
C ZERO ARRAYS N 15 

693 DO 1 J=l t 20 N 16 

694 DATA (J)=0. N 17 

695 1 CFS(J=0. N 18 

696 ID1=1 M 19 
C FIND RATING CURVE WITH SMALLEST MAXIMUM FLOW RATE N 20 

697 2 QMIN=QI20,IOU M 21 

698 MIN=ID1 N 22 

699 GO TO 4 N 23 

700 3 ID1ID1H N 24 

701 IF (QMIN-Q{20 t IDU) 4,4,2 N 25 

702 4 IF tlOl-NOVS) 3,5,5 N 26 

703 5 1=1 M 27 



704 oo 6 J = 2,20 Y l "~ LOHEST RATING CURVE * 28 

705 

706 6 

C 

707 
708 
709 
710 

711 7 

7L2 17 

7L3 
714 
71 5 8 

717 

718 9 

719 

720 

C 

721 
722 
723 
724 
725 

727 10 

C 

726 
729 11 



, z 

C 
732 L3 



733 U 

734 15 



735 16 

736 

737 



"" 



739 
740 
741 
712 
7<i3 
744 

r c 

T46 I 



SCFSI I )=J( J,MIN) 


N 


29 


1=1+1 


N 


30 


COMPUT END AREA AND DEPTH 


N 


31 


(30 9 ID1=1,NOVS 


N 


32 


DO 9 J=l,19 


N 


33 


00 7 I =2,20 


N 


34 


IF <Q< I,ID1)-SCFS(JM 7,17,8 


N 


35 


CONTINUE 


N 


36 


DATA ( J)=A( I, IOl)+DATA{J) 


N 


37 


CFS<JI=OEEP(I .IDU + CFSU) 


N 


38 


GO TO 9 


N 


39 


DATA ( J)=A( I- 1, II) U *XY*( A( I 101) I f T I in? 1 '* 

CFS<J)=OEEPU-l,IDl)tXY*(DEEPII,ID 1 ncpoi T i ?I?[ ^ 
CONTINUE <"ttf 11,101 ) -DEEP I 1-1, ID1 M fCFS( J) 


N 
N 


:i 


XNCVS=NQVS 


N 


44 


WRITE (6,13) REACH 


N 


45 


COMPUTE TRAVEL TIME 


N 


46 


DO 10 1=1,19 


N 


47 


AVAREA=DATA{ I J/XNQVS 


N 


48 


OPU) = CFS( IJ/XNQVS 


N 


49 


S=AVAREA*XLD36 


N 


50 


c* n=s/scFS(i j 


N 


51 


WR ITE (6,141 DP(I|,SCFSU).CU> 


N 


52 


CONTINUE 


N 


53 


PUNCH CODE 


N 


54 


IF (NPU) 12,12,11 


M 


55 


WRITF, (7,15> ID, REACH, XL, SLOPE 


N 


56 


RETURN 


N 


59 


FORMAT<1HO,T46, 'TRAVEL TIME TABLE-/T54 iRFArwi ^ *,,-,,, 
156, 'FLOW T65r'TRAVEL'/T46,'DFPTH' T56 IRATP" rl? " //T * 6 r W AT ER ' ' T 


N 


61 


2' ,T56, 'CFS' ,T66,'HRS') * ? ? 1 &6t ' TIME /T46, FEET 


N 


62 


FORMAT UOX,F10.2,FIQ.O,F10.4) 


n 


63 


FORMAT( 'STORE TRAVEL T IVf ' , \? 1 , ' in- . n 770 or.sru ., n 
1' LENGTH^, F9.0,' FT ' /T2 1 , ' SLGPE = ! t F all I J??FT? N0= f F5. 1 t T44, 


N 
M 


65 


21,'DEPTH(FT)*T35i*FLOW(rP^M T/O iriur-* * ' ' /T2 


N 


66 


FORMAT IT21 t F7.2,Fl5.0,F15.3) ' ' ' TIMECHRS)r| 


N 


67 


END 


N 


68 




N 


69- 


SUBROUTINE ROUTE 





i 










1 


ID=DATA{1) 





8 


NHODATA(2 





9 


IDH=DATA(3) 





10 


OT{IO)=DATA(4) 





U 


DA(ID)=OA( [DHJ 





12 


M=IENO( I OH) 


a 


13 


IF ID AND IOH ARE EQUALt ADD I TO I DH 


a 


14 


IF (IDIDH) 3,1,3 





l& 


IOh=IDH+l 





16 




a 


17 
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747 DO 2 1=1, M 

7*8 2 OCFS( I iIDH)OCFS(IiIDH-l) 

749 DT( IDH-DT{ IOH-1) 

750 REAM IDH = PEAK UDH-U 

751 3 NERRT=0 

752 PEAK(ID)=1. 

753 R0=0. 

754 N=19 

755 OCFSU,IO)=0. 

756 5=0. 

757 T 1=0(1) 

758 J=l 

759 GUES-l. 

760 CFStl)=0. 

C IF ROUTING INTERVAL IS NOT EQUAL TO TIME INCREMENT OF INFLOW 

C HYDROGRAPH, I INTERPOLATE 

761 IF (DTUD)-OT( IDHM 8,15,4 

762 4 TID=DTUD) 

763 TIDH=0. 

764 DD 7 1=2, M 

765 TIDH=T IDH+DTt IDH) 

766 IF (TID-TIOHI 6,5,7 

767 5 J=J+L 

768 CFS{ Jl=OCFS( I , IDH) 

769 TIC-TID+OT( ID) 

770 GG TO 7 

771 6 J=J+l 

772 CF5( J)-OCFSU-ltIDHI + < t TID-T IOH + OT ( IDHH /DT( I DH) I * ( OCFS ( I , I DH)~DC I 
1S( I-li IDHJ ) 

773 TIC=TID+DT(IDJ 

774 7 CONTINUE 

775 GG TO 13 

776 8 TIDH=0. 

777 TID=OTtID 

778 DO 12 1=2, M 

779 TIDH=TIDm-DT( IDH) 

780 9 IF (TIDH-TID) 12,10,11 

781 10 J*J+1 

782 CFS( J)=OCFS( I , IDH) 

783 TI D=TID+DT( ID) 

784 IF (J-3001 12,13,13 
7B5 II J=J+1 

786 CFS(J)=OCFS( 1-1, IDH)*f ( TI D-T IDH+DT ( I DH) ) /DT I I DH 1 >* t OCFS (I , I DH ) -OC \ 
1SI 1-1, IDH) ) 

787 TIC=TID*DT(IDI 

788 IF (J-300) 9,13,13 

789 12 CONTINUE 

790 13 IEND(IDH=J 

791 OT( IDH)=UT( ID) 

792 M=J 

793 DO 1.4 1 = 2, M 

794 14 OCFSI I , IDH)=CFStI ) 

C IF INFLOW IS ZERO, SO IS OUTFLOW 

795 15 DO 16 L=2,H 

796 IF (OCFStL, IDH) ) 16,16,49 

797 16 OCFSIL ,ID)=0. 
C ROUTE 

798 49 DATA (L-l)=0. 

799 DO 42 I=L,300 
000 IF (I-M) 18,18,17 



801 17 OCFS(I,IDH)=DCFS( I-l.tDHI* 9 

18 



804 J=1 RO 



so, - s92 

" I -' SM S 83 



R4 

808 GO TO 25 n RC 

809 20 DO 21 J = 2,N D Q 6 

810 IF <"CFS([,IDHI-SCFS(JII 24,23.21 87 
flll 21 CONTINUE t^3,21 r 

812 If= (NHRRT) 22,22,36 G 89 

813 22 WRITE (6,46) n 

814 NERRT=1 I, 

815 GO TO 36 n 97 
S16 23 DI2=DP(J) <^ 

817 TI2=C(JI g4 

818 GO TO 25 9 g 

24 





823 J=l 



TIHE OF OUTFLn S iSi 

825 26 D02=(GUES/SCFS(l))*DP(l 103 

826 T02 = Cm , 04 

827 GO TO 31 10 g 

828 27 DO 2fl J=2,N 106 

829 IF (GOES-SCFS(J)J 30,29,28 107 

830 28 CONTINUE 108 

831 J=N lrt g 

832 29 DQ2 = DP(J) ^Q 
33 T02=C( J) 111 
S34 GO TO 31 112 

835 30 RA no=<GUF.S-SCFS(J-lM/(sCFS(J)-SCFSli-in 

836 =- 



_ 

837 TOa=C(J-l)+RATIO*IC(J)-CtJ-l)| 115 
C FIND WATER SURFACE SLOPE 116 

838 31 DDD-DISTIIDI/(DIST(ID)+C[2-D02) U7 
639 IF (DDD-.Ol) 32,32,33 118 
840 32 GUES=OCFS( I-l.IDHJ 119 
8*1 GO TO 35 120 
842 33 T2=.5*(TI2tTU2) 121 
8 ^ 3 T2=T2*SQRT(DOO) 122 
8** T = TUT2 123 

C COMPUTE ROUTING CQEFFIUFNT 124 



816 Q2=COF.F*SIA 126 

847 TRY1=GUES G 127 

848 RATIO=0?/(GUES+IE-20 128 
3^9 OIFF=AHS(1. -RATIO) n 129 

C TEST FOR CONVERGENCE 130 

850 IF (DIFF-.001) 37,37,34 131 

851 34 GUES=D2 132 
^52 35 CONTINUE 133 
853 UCFS<I,II)*DATA(I-1|*S[A 134 
a54 DATA ( I)=[JATA( 1-1) 135 
8 " WRITE (6,47) I t nCFS(I,!D) n & 

137 
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856 GO TO 38 



138 



857 36 QCFSUiIUI-DATAU-mSIA l39 

858 DATA ( I )=OATA ( I -1 ) l/f 

859 GO TO 38 141 

860 37 OCFSU, 101=02 n l42 

861 DATA I I )=CQEF 1 *3 
C COMPUTE NEW STORAGE 14* 

862 38 S=SIA-OCFSU, ID) 1*5 

863 Tl=T2 D U6 

864 RO=RQtQCFS( I, ID) a l^ 7 

865 IF (OCFStI tIDJ-OCFSU-lf IOJ 39,40,40 148 

866 39 IF lOCFSt I,ID)-1.) 43,43,42 149 

867 40 IF (OCFSU ,IO)-PEAK(IDM 42,42,41 150 

868 41 PEAK( IO)=t)CFS( I,ID) 151 

869 42 CONTINUE 15 2 

870 1=300 153 

871 43 IEND( !D) = I n 15/t 

872 ROIN=1RO*DT(IO) )/(DA( 101*645.333) 155 
C PUNCH CODE 156 

873 IF (NPU) 45,45,44 D 157 

874 44 WRITE (7,48) ID ,NHD , IQH ,DT( I D) 158 

875 45 RETURN 159 
C 160 

876 46 FORMAT1 1HO, 'TRAVEL TIME TABLE EXCEEDED') 161 

877 47 FORMATIT10, 'PROBLEM FAILED TO CONVERGE AFTER10 ITERATIONS. CONVERG 162 

1ENCE WAS FORCED. '/T20, 'OUTFLOW NUMBEK = ',14, 'RATE =',F10.2) 163 

878 48 FORMAT! ' ROUTE ,T21 , ' 10 = ' , 1 1 f T29, ' HYD N0= ' ? I 3 ,T45, ' INFLOW ID=',I 164 

H,T65 f < DT= 1 ,F8.6,'HRS'l 165 

879 END 166- 

880 SUBROUTINE RESVQ P 1 
C THIS SUBROUTINE ROUTES A HYDROGRAPH THROUGH A RESERVOIR WITH THE P 2 
C STORAGE-INDICATION METHOD. P 3 

881 COMMON CFS(300) , OCF S ( 300,6 ) , I END (6 ) ,OATA t 3 10 , DA< 6 ) r DP ( 20 ) ,NPU,NHD P 4 
1, SCFS( 20) ,C( 20) ,A(20, 6) ,Q(20,6) ,RAIN( 200) , DEEP t 20,6) ,NER , MAXNG, NCO P 5 
2MM,ICC,NCODE,DIST(6),SEGN(6),CTBLE(50,11) , ITBLE ( 50 , 2 ) , ZALFA ( 20 ) , DT P 6 
31 6),TIME,PEAK(6)tROIN,ISG(6) P 7 

882 ID=DATA(1) P 8 

883 NHO=DATA(2) P 9 

884 IOH=DATA(3) P 10 

885 NERES=0 P 11 

886 DTUD) = DT( IDH) P 12 

887 RQO. P 13 

888 DA(ID)=DAI IDH) P 14 

889 PEAK(ID)=1. P 15 

890 J=l P 16 

891 1-4 P 17 
C REGAINING DATA ARE FLOW AND STORAGE VALUES P 18 

892 SCFS(J)=DATA( I) P 19 

893 STCREl=OATA(m)*12.1 P 20 

894 STORE=STORE1 P 21 
C COMPUTE STORAGE COEFFICIENT ARRAY C P 22 

895 1 C( J)=(SCFS(JI/2.)+(STORE/DT(ID)) P 23 

896 1 = 14-2 P 24 

897 J=J+1 P 25 

898 IF (J-20) 2,2,3 P 26 

899 2 SCFS(J)=DATA(I) P 27 

900 STORE=DATA(H-1)*12. 1 P 28 

901 IF (SCFS( J)-.001) 3,3,1 P 29 

902 3 N-J-1 P 30 
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903 OCFSU, 101=0. 

904 s=STrwci/nnici p 31 

C ROUTE p 3 2 

905 DO IS 1=2, 150 P 3-5 

906 IF f I-IEND< IOH) ) 5,5,4 P 3* 

907 4 OCFSU, !f)H)=0,0 P 35 

908 "5 AVIN=[nCFSU,IOHUQCFS<I-l t i DH M/ ? P 36 

909 SIA=S+AVIN *t 1UHM/2. p ^ 

C DETERMINE PROPER C P 38 

910 DC) 6 J = 1,N P 39 
9H IF ISIA-UJM 10,9,6 P 40 

912 6 CONTINUE P 41 

913 II- (MF-HCS) 7,7,8 P 4? 

914 7 URITE if,, 19) p 43 

915 NRRf:S=l p 44 

916 RESC=SChS(N)/C<N P 45 
C COMPUT OUTFLOW P 46 

917 OCFSII, ID)=RGSC*SIA P 47 

918 GO TO 11 p 48 

919 9 OCFSU, ir=SCFS(J) P 49 

920 GO TO 11 P 50 

921 10 I'ICFSl I, 10)=SCF5( J-lH-miA-ri 1-1 i i/ /r/ .. f , p 51 

11)) 5[A C(J 1)l/(CtJ '-C(J-l)), MS CFS(J)-SCFS(J- P 52 

C DETERMINE NEW STORAGE P 53 

922 11 S = SIA-fJC(-S( I r IU) P 5.4 

923 KO-Rl)K)CI : SU, IDI P 55 

924 IF (OC.FSf I,I[))-QCI-SU-l t lD)) 12,13.M P 36 

925 12 IH ((.)CFSU,IDI-M 16,16,15 1 "' IJ ' 13 P 57 

926 13 IF lUCFStMOI-PFAKHDM 15,15.14 P *>* 

927 14 PFAK(II,l) = nCFS(I ,![)) ^tit p 5g 
92B 15 CONTINUE P 60 

929 I-lliO P 61 

930 16 IFNI)(IDI = I P 62 

931 K[)1N"RIJ*TU|))/(I)A< I [) ) *645. 333 1 p ^3 
C PUNCH CODl- P 64 

932 II-' (NPU) 1H,1B,17 P 65 

933 17 WK ITl: (7,201 I1),NNIJ, IDH p 66 

934 I I?.*Ni-3 P 67 

935 HIUTI: (7,211 UJATAU , [a5, 1 I I P 68 

936 IB RF.TURN P 69 
C P 70 

937 I'J FORMAT ( 1110 , 3 3MSTORAGK-0 1 SCHARGE TABLE EXCEEDFO 1 P 71 

938 2 ^r^.u"^^ p 7 2 3 

ZAGtitAC NT)M /T21 ,' OUTFLOW (CFS) ' ? T37, STOR P 74 

939 21 FOUMAT ( T2 1 ,f : 10, 1 ,F 1 3. II P 75 

940 END p 75 

P 76- 

941 SUIWOUTINI: UKROR 

C 



o 5 

* 

It)lOATA(ll 7 



945 ssro. 9 

*& WRirii (6,15) J? 

947 ji Q 11 

12 



73 



C IF TIME INCREMENTS NOT EQUAL, INTERPOLATE 13 

948 IF (DTUDU-DTU02M 1,8,2 14 

949 1 L-1DI O }? 

950 K=ID2 n }? 

951 GO TO 3 Q -|5 

952 2 L=ID2 n fa 

953 K=I01 ,q 

954 3 M=IEND(L) n 

955 TIO=DT(K) n ?V 

956 TIDH=0. yy 

957 DO 6 1=2, M \\ 

958 TIDH=TIDH*OTU) n yl 

959 IF (TID-TIOHI 5,4,6 n ?K 

960 4 j=j+-i ;; " 

961 CFSI J)=OCFSU 



962 TID=TID+OT(K) " X Q II 

963 GO TO 6 i 



5 J = J+1 

965 CFS(jrOCFSU-l,L)t'((TID-TIDH*DT(LI )/DT(L) )*(OCFS( I , D-OCFS ( 1-1, L) Q 31 

966 TIC-TID+OTU) S ?! 

967 6 CONTINUE n ff 

968 IENOUUJ 

969 DT(L)=DT(K) ,' 

970 DO 7 I-2.J O ^? 

971 7 OCFSU,U = CFSU) n fg 

972 8 IF (IENOI 1 01 )-I END [ ID2 ) ) 9,9,10 Q ^q 

973 9 M=IEND(IDi) n In 

974 GO TO 11 " 7" 

975 10 M=IEND(ID2J JJ 71 

976 11 T2=TIME 
C DETERMINE ERROR n ll 

977 DO 12 1=1, H J3 

978 ERR=QCFSU, ID 1)-OCFS( 1,1021 n 46 

979 ^RITE (6,16) T2,OCFSU,IDn,OCFS(I,I02),ERR 47 

980 T2=T2+DT(ID1) J J' 
C SUM OF SQUARES OF ERROR ?q 

981 12 SSE=SSE+ERR*ERR n n 

982 XM=M " 
C ERROR VARIANCE n l-> 

983 EVAR=SSE/XM J ff 
C ERROR STANDARD DEVIATION O .54 

984 ESDEV^SQRT(EVAR) ll 

985 WRITE (6,17) ESOEV Q II 
C PERCENT ERROR FOR PEAK FLOWS n c 7 

986 ERPK=ABS(PEAK(ID1I-PEAK(ID2I) 5 50 

987 PCTER=(EKPK/PEAK(ID1))*100. n cq 

988 WRITE (6,18) PCTER n ^n 
C PUNCH CODE Q ?, 

989 IF (NPUI 14,14,13 I- 

990 13 WRITE (7,19) ID1,ID2 O h l 

991 U RETURN J ^ 

992 15 FORMAT! 1HD.T33, 'TIME' ,T55, 'FLOW l,T76 f i F LQW 2 ' , T95 , ' ERROR' /T34, Q 66 

l'HRS',T57,'CFS',T78,'CFS',T97,'CFSM n 67 

993 16 FORMAT (20X , F20 . 3 , 3F20.0) 68 

994 17 FORMATI1HO,T10, 'ERROR STANDARD DEVIATION = ',F10.3 Q 69 

995 18 FORMAT(T10,'PEAK DISCHARGE ERROR = ',F7.2,' PERCENT'///) Q 70 

996 19 FORMAT! 'ERROR ANAL YS I S ' , T21 , ' ID I=',I1,T29,MD II=',I1( Q 71 

997 END 
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SUBROUTINE SEDT 

THIS SUBROUTINE COMPUTES THE SEDIMENT YIELD FOR A FLOOD 

COMMQN CFS(300) ,OCFS(300,6) , 1 END! 6 ) ,DATA ( 3 10) t OA( 6 ) , DPI 20 ) ,NPU ,NHD 
1, SCFSUO),C(20) ,A(20,6) ,Q (20 , 6 I , RAI N I 200 ) , DEEP ( 20 , 6) ,NR , MAXNO.NCQ 
2MM 1 ICCiNCODE,DIST(6ltSEGN(6lfCTBLE(50U)ITBLE<50,2),ZALFA120),DT 
31 6),TIME,PEAK(6, RQINi ISG(6) 

IO=DATA( II 

SOIL=DATA( 2) 

CROP=OATA{3) 

CP=DATA(4) 

SL=DATA(5) 

COMPUTE SEDIMENT YIELD 

X=ROIN*DA( rot*53.333#PEAK(IDI 

SED=95.*X**.56*SOIL*CRGP*CP*SL 

WRITE (6,3) SED 

PUNCH CODE 

IF INPUI 2,2,1 

WRITE (7,4) ID, SOILtCROP,CP,SL 

RETURN 



FORMAT HOX 
FORMAT ( 
1=' ,F5.3 
END 



'SEDIMENT YIELD = 



F10.lt ' TONS') 



Xr 'SEDIMENT YIELD = , F.t 

'SEDIMENT Y IELD * , T21 , ID = ' , II ,T29, ' SOIL= ' ,F5.3,T42 , ' CROP 
P = ',F5.3,T70,'LS*',F5.3) 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24- 



//SDATA 
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ZALFA = 1234567890 #. , - 



COMMAND TABLE 



START L 2 

STORE HYD 2310 

RECALL HYD 3310 

COMPUTE HYD 4310 

PRINT HYO 5 2 

PUNCH HYD 6 1 

PLOT HYD 7 2 

AOD HYD fl ^ 

STORE RATING CURVE 9100 

COMPUTE RATING CURVE 10310 

STORE TRAVEL TIME 11100 

COMPUTE TRAVEL TIME 12 5 

ROUTE 13 4 

ROUTE RESERVOIR 14100 

ERROR ANALYSIS 15 2 

SEDIMENT YIELD 16 5 

FINISH 17 
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